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Abstract 


A measurement is presented of the tt inclusive production cross-section in pp collisions at 
a center-of-mass energy of = 8 TeV using data collected by the ATLAS detector at the 
CERN Large Hadron Collider. The measurement was performed in the lepton-i-jets final state 
using a data set corresponding to an integrated luminosity of 20.3 fb”*. The cross-section 
was obtained using a likelihood discriminant fit and ^-jet identification was used to improve 
the signal-to-background ratio. The inclusive ft production cross-section was measured to be 
260 + l(stat.)t 23 (sysb) + 8(lumi.) + 4(beam) pb assuming a top-quark mass of 172.5 GeV, in 
good agreement with the theoretical prediction of 253;|;[ j pb. The tt —> {e, p) -i-jets production 
cross-section in the fiducial region determined by the detector acceptance is also reported. 
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1 Introduction 


The measurement of the top-quark pair production cross-section is an important part of the physics pro¬ 
gram of the Large Hadron Collider (LHC) [1]. As the top-quark is the heaviest known fermion, it has the 
largest coupling to the recently discovered Higgs boson and plays a special role in many theories beyond 
the Standard Model (SM). New physics may result in additional top-quark decay channels or new mecha¬ 
nisms of W production that can cause the measured cross-section to deviate from the SM prediction. Also, 
tt production is the dominant background to many searches for new physics. 

The inclusive pp ^ tt cross-section cr,,- at a center-of-mass energy of = 8 TeV has been calculated 
to be 253^}j pb for a top-quark mass of mtop = 1V2.5 GeV. These calculations were performed at next- 
to-next-to-leading order (NNLO) in QCD including resummation of next-to-next-to-leading logarithmic 
(NNLL) soft gluon terms using the program top-r-i- 2.0 [2-7]. In these calculations, the renormalization 
scale jUr and factorization scale were both set equal to witop. The NNLO-i-NNLL value was found to be 
about 3% larger than the NNLO-only calculation as implemented in Hathor 1.5 [8]. As a result of recent 
progress in theoretical calculations, the tt production cross-section has been evaluated with uncertainties 
below 6%. In calculating these uncertainties, the scale uncertainty, evaluated by independent variations of 
/tr and jUp by factors of 1 /2 and 2, was combined in quadrature with uncertainties on the coupling strength 
of the strong interaction, Q's^ and the parton distribution functions (PDFs). The latter were calculated 
using the PDF4LHC prescription [9] with MSTW2008 NNLO (at the 68% confidence level) [10,11], 
CTIO NNLO [12,13] and NNPDF2.3 5f FFN [14] PDF sets. 

In the SM, the top-quark decays into a VT-boson and a f7-quark with a branching ratio close to 100%. 
The present analysis aims to measure the tt production cross-section in the lepton-i-jets final sfafe, where 
one of fhe IT-bosons decays info an elecfron or a muon (collecfively called a lepton and denofed i) and 
a corresponding neufrino (including fhe W tVt eVrVe and W rvr ^ decays), and fhe 

ofher VF-boson decays hadronically. The final sfafe is characferized by fhe presence of a highly energefic 
isolated lepfon, large missing fransverse momenfum (F™®®) due fo fhe neulrino(s) escaping defecfion, and 
four jefs due fo fhe fwo 6-quarks from fhe fop-quark decays and fhe fwo quarks from fhe hadronic W 
decay. The selecfed evenfs are required fo have af leasf fhree jefs, allowing one Jef fo be undetected due fo 
limifed defector accepfance and jef reconsfrucfion inefficiency. Af leasf one of fhe jefs is furfher required 
fo be idenfified as a 6-jef. 

The inclusive fop-quark pair production cross-section in pp collisions af 8 TeV has been previously mea¬ 
sured in fhe dilepfon channel by bofh fhe ATLAS and CMS Collaborations [15, 16]. The resulfs are 
found fo be in good agreemenf wifh fheorefical predicfions. This paper presenfs fhe measuremenf of fhe tt 
producfion cross-secfion of fhe ATLAS Collaborafion in fhe lepfon-i-jels channel, which provides a cross¬ 
check of fhe dilepfon measuremenf fhaf has differenf background conditions and sysfemafic uncerfainfies, 
and confribufes fo fhe combined tt producfion cross-secfion resulf. In addition, fhe measuremenf in fhe 
lepfon-i-jefs channel is imporfanf for probing fhe presence of new physics fhaf changes fhe fop-quark decay 
branching fractions, e.g. producfion of charged Higgs bosons via t —> H^b decays. 

The measuremenf was performed using fhe full 8 TeV dafa sef (20.3 fb“^). In addition to fhe fofal tt 
producfion cross-secfion, a fiducial cross-secfion was measured, defined using physics objecfs consfrucfed 
of sfable particles fo approximafe fhe tt 6y-i-jels defecfor accepfance. 

The majorify of fhe evenfs in fhe selecfed sample originafe from tt producfion. However, evenfs from 
ofher SM processes (IT/Z-i-jefs, single fop-quark, diboson, and mulfijef producfion) are also presenf. To 
defermine fhe signal fracfion, a discriminating variable (likelihood discriminanf, LHD) was consfrucfed 
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based on the kinematic properties of Monte Carlo (MC) simulated signal and background events. A 
weighted sum of LHD distributions (“templates”) for the signal and for the background was fitted to the 
LHD distribution in data, and the resulting number of signal events was converted to the tt inclusive 
production cross-section using the signal-reconstruction efficiency (determined from MC simulation), 
known branching ratio for the lepton-i-jets final state, and the total integrated luminosity. 

The paper is organized as follows. The ATLAS detector is briefly described in Sec. 2, followed by 
descriptions of the data and MC samples used in the analysis (Sec. 3). The event selection and the 
backgrounds are outlined in Sec. 4. The measurement procedure is presented in Sec. 5 for the inclusive 
case, and in Sec. 6 for the fiducial case. The estimation of the systematic uncertainties is presented in 
Sec. 7. Finally, the results and conclusions are given in Sec. 8 . 


2 The ATLAS detector 

The ATLAS detector [17] at the LHC covers nearly the entire solid angle around the collision point. It 
consists of an inner tracking detector surrounded by a thin superconducting solenoid, electromagnetic 
and hadronic calorimeters, and a muon spectrometer incorporating three large superconducting toroid 
magnets. The inner-detector system (ID) is immersed in a 2 T axial magnetic field and provides charged 
particle tracking in the range I 77 I < 2.5. 

The high-granularity silicon pixel detector is closest to the interaction region and typically provides three 
measurements per track, the first hit being normally in the innermost layer. It is followed by the silicon 
microstrip tracker which has four layers in the barrel region. These silicon detectors are complemented by 
the transition radiation tracker, which enables radially extended track reconstruction up to I 77 I = 2.0. The 
transition radiation tracker also provides electron identification information based on the fraction of hits 
(typically 30 in total) above a higher energy deposit threshold corresponding to transition radiation. 

The calorimeter system covers the pseudorapidity range \r]\ < 4.9. Within the region \t]\ < 3.2, electro¬ 
magnetic calorimetry is provided by barrel and endcap high-granularity lead/liquid-argon (LAr) electro¬ 
magnetic calorimeters, with an additional thin LAr presampler covering \r]\ < 1.8, to correct for energy 
loss in material upstream of the calorimeters. Hadronic calorimetry is provided by the steel/scintillating- 
tile calorimeter, segmented into three barrel structures within \r]\ < 1.7, and two copper/LAr hadronic 
endcap calorimeters. The solid angle coverage is completed with forward copper/LAr and tungsten/LAr 
calorimeter modules optimised for electromagnetic and hadronic measurements respectively. 

The muon spectrometer (MS) comprises separate trigger and high-precision tracking chambers measuring 
the deflection of muons in a magnetic field generated by superconducting air-core toroids. The precision 
chamber system covers the region \t]\ < 2.7 with three layers of monitored drift tubes, complemented 
by cathode strip chambers in the forward region, where the background is highest. The muon trigger 
system covers the range \r]\ < 2.4 with resistive plate chambers in the barrel, and thin gap chambers in the 
endcap regions. A three-level trigger system is used to select interesting events [18]. The Level-1 trigger 
is implemented in hardware and uses a subset of detector information to reduce the event rate to a design 
value of at most 75 kHz. This is followed by two software-based trigger levels which together reduce the 
event rate to about 200 Hz. 
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3 Data and simulated samples 


The data were collected during the 2012 LHC running period at a center-of-mass energy of 8 TeV. After 
applying data-quality selection criteria, the data set used in the analysis corresponds to an integrated 
luminosity of 20.3 fb“^ 

Simulated tt events were generated with Powheg [19] interfaced to Pythia v6.426 [20] for the fragmenta¬ 
tion and hadronization, with the next-to-leading order (NLO) CTIO PDF set. The Perugia2011C underly¬ 
ing event tune [21] with the CTEQ6L1 PDF set [22] was used. An alternative ttMC sample was produced 
with MC@NLO v4.01 [23] using the CTIO PDF set. In this sample, the parton shower and the under¬ 
lying event simulations were performed with Her wig v6.520 [24] and Jimmy v4.31 [25], respectively, 
using the AUET2 tune [26]. To estimate the model dependence of the parton shower and fragmentation 
modeling, the signal sample generated with Powheg interfaced to Pythia was compared to that generated 
with Powheg interfaced to Herwig-eJimmy. Additional tt samples simulated with AcerMC [27] interfaced 
to Pythia based on the CTEQ6E1 PDE set were used to study the systematic uncertainties arising from 
initial- and final-state radiation. All tt samples were produced with a top-quark mass of 172.5 GeV and 
normalized to the NNEO-i-NNEE cross-section quoted in Sec. 1. 

The dominant background to tt production is vector-boson production in association with jets, F-i-jets 
{V = W,Zly*). Samples of events were generated using Alpgen [28] interfaced to Pythia v6.426 based 
on the CTEQ6E1 PDE set and the Perugia2011C tune. The MEM parton and jet matching procedure [29] 
was applied inclusively for F-E5-light-partons (2^7) production and exclusively for lower multiplicity 
samples. In addition to F-i-light partons, the production of vector-bosons with additional heavy-flavor 
partons (V + c-i-jets, V + cc-i-jets, V + bb+iets) was also simulated. Inclusive F-i-jets samples were formed 
by combining the samples of light and heavy quarks according to their respective cross-sections. An 
alternative VT-i-jets sample simulated with Sherpa [30] with massive fi/c-quarks was produced to evalu¬ 
ate systematic uncertainties on VF-i-jets modeling. The F-f-jets processes were normalized to the NNEO 
calculations [31,32]. 

The background contribution from single-top production, including t- and ^-channel contributions and 
Wt production was simulated with Powheg interfaced to Pythia using the CTIO PDE set. Einally, diboson 
production (WW, VFZ, ZZ) was simulated with Herwig using the CTEQ6E1 PDE set. All samples were 
scaled according to their theoretical production cross-sections: at approximate NNEO for the single top t- 
channel [33], and at NEO-i-NNEE for the single top ^-channel [34] and ITt-production [35]. The diboson 
processes were normalized to the NEO [36] predictions. 

Most of the tt samples and all the background samples were processed through the full ATEAS detector 
simulation [37] based on GEANT4 [38]. A few tt samples used for the evaluation of certain systematic 
effects (modeling of initial- and final-state radiation and parton showers - see Sec. 7 for details) were 
produced using the ATEAS fast simulation that employs parameterized showers in the calorimeters [39]. 
All MC samples were reconstructed using the same analysis chain as used for the data. Correction factors 
were applied in order to better reproduce the trigger, lepton reconstruction efficiency, and b-]et idenfifica- 
fion efficiency observed in fhe dafa. The simulations also included fhe effecl of multiple pp collisions per 
bunch crossing (pileup). 
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4 Event selection 


The events were required to pass a logical OR of isolated and nonisolated single-lepton {e or /i) trigger 
conditions with transverse momentum px thresholds of 24 GeV for isolated triggers and 60 (36) GeV 
for nonisolated single-e (single-p) triggers. All events were also required to have a reconstructed hard 
collision primary vertex (the main PV) built of at least five particle tracks with px > 0.4 GeV. 

The reconstructed objects used in the analysis include electrons, muons, jets, and missing transverse 
momentum. Electron candidates were reconstructed from an electromagnetic energy deposit matched 
to a track in the inner detector [40]. They were required to have transverse energy Ex > 40 GeV and 
pseudorapidity of the cluster^ \r]\ < 2.47, excluding the barrel-endcap transition region 1.37 < \r]\ < 1.52. 
Electron candidates were also required to originate from less than 2 mm along the z-axis from the main 
PV and satisfy isolation criteria. The latter involve a combination of calorimeter isolation (a requirement 
on a sum of energies of calorimeter cells within a cone of size A/? = 0.2 around the electron direction) 
and track isolation (a requirement on a scalar sum of track px within a cone of size AT? = 0.3, in each 
case excluding the contribution from the electron itself). Both requirements were chosen to separately 
result in a 90% electron reconstruction efficiency for prompt electrons from Z ^ ee decays. Muon 
candidates were reconstructed using combined information from the muon spectrometer and the inner 
tracking detectors [41]. They were required to have px > 40 GeV and Ipl < 2.5, and, like electrons, to 
originate from within 2 mm along the z-axis of the main PV. The muon isolation was defined in ferms of 
fhe rafio of fhe scalar sum of fhe frack px in a cone of variable radius A/? = 10 GeV/pf^ around fhe muon 
direction (excluding fhe muon frack if self) fo fhe px of fhe muon (p^). This rafio was required fo be less 
fhan 0.05, corresponding fo a 97% selecfion efficiency for prompf muons from Z —> pp decays. 

Jefs were reconsfrucfed using fhe anfi-kj algorifhm [42,43] wifh radius paramefer R - 0.4. The measured 
jef energy was correcfed for inhomogeneifies and for fhe noncompensafing nafure of fhe calorimefer 
fhrough jef px- and p-dependenf factors derived from MC simulation and validated wifh dafa. Jefs were 
calibrated using a combination of in-sifu fechniques based on fhe fransverse momenfum balance befween 
a jef and a reference objecf [44]. The jefs were required fo have correcfed px > 25 GeV and |p| < 2.5. 
Jefs from pileup were suppressed by requiring fhe absolute value of fhe jef vertex fracfion^ for jefs wifh 
Px < 50 GeV and Ipl < 2.4 fo be above 0.5. 

Jefs were fj-fagged (idenfified as originafing from fj-quarks) using fhe mulfivariafe-based algorifhm algo¬ 
rifhm (MVl) [45]. This is a neural nefwork-based algorifhm fhaf makes use of frack impacf parameters 
and reconsfrucfed secondary verfices. Jefs were idenfified as f7-jefs by requiring fhe MV 1 oufpuf fo be 
above a cerfain fhreshold value. This value was chosen such fhaf fhe overall lagging efficiency for b- 
jels originafing from lop-quark decays in MC tt evenls is 70%. Tagging scale faclors were applied fo 
correcl for fhe difference in fhe lagging efficiency befween MC simulation and dafa, including fhe inef¬ 
ficiency scale faclors applied when a jef was nol lagged. The scale faclors were derived as explained in 
Refs. [45-47]. 

To avoid double-counfing objecfs in an even! and fo suppress leplons from heavy-flavor decays, fhe fol¬ 
lowing procedure fo remove overlaps was applied to fhe reconsfrucfed objecfs: (1) removal of jefs malched 
wifhin A/? = 0.2 of elecfrons wifh Ex > 25 GeV; (2) removal of eleclrons malched wifhin AE = 0.4 of 


* Unlike the (regular) pseudorapidity determined using the object direction, the pseudorapidity of the cluster is defined using 
the position of the reconstructed cluster in the calorimeter with respect to the geometric center of the detector. 

^ The Jet vertex fraction is defined using the tracks matched to a jet as the ratio of the scalar sum of the transverse momenta of 
tracks from the main PV to that of all tracks. A jet without any track matched is assigned a jet vertex fraction value of -1. 
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jets with pj > 25 GeV; (3) rejection of events where the selected electron shares an inner detector track 
with a selected muon; and (4) removal of muons matched within AR = 0.4 of jets with pj > 25 GeV. 
After these steps, the events were required to have exactly one lepton (an electron or a muon) selected 
as above and matched to the trigger object that caused the event to be selected, and at least three jets, 
of which at least one was (7-tagged. Events with a second lepton with pj > 25 GeV were discarded.^ 
Finally, the events had to exceed minimum values for the W transverse mass mxClV) and the magnitude of 
the missing transverse momentum The value of in the tt —> ^v+jets process is relatively large 

due to neutrino(s) from W decays. It also includes energy losses due to detector inefficiencies as well 
as energy fluctuations in jet measurements. The missing transverse momentum value is calculated as the 
negative of the vector sum over the energies of all clusters in the calorimeters, and is refined by the 
application of the object-level corrections for the contributions arising from identified electrons, muons, 

and jets. The VT-boson transverse mass mjiW) is defined as mT'{W) - yj2p^pj{\ - cos{ip^ - (p'')), where 

Pj, and p^, refer to the transverse component and ip value of the lepton momentum and the missing 
momentum, respectively. The events in this analysis were selected by requiring E™'*® > 30 GeV and 
nniW) > 30 GeV. 

After applying the selections listed above, the signal efficiency for tt events with at least one top-quark 
decaying leptonically was found to be 4.4% in the e-i-jets channel and 5.5% in the p+jets channel. 

While the majority of events in the selected sample include ‘real’ leptons (prompt electrons and muons 
from vector-boson decays), a small fraction of events is due to leptons referred to as ‘fakes’ (which include 
nonprompt leptons and artifacts of the reconstruction). There are several sources of fake leptons. Fake 
muons predominantly come from semileptonic decays of heavy-flavor hadrons. In the e-i-jets channel, in 
addition to this source, there are electrons originating from conversions as well as photons and hadrons 
misidentified as electrons. The signal sample is composed of events with either real or fake leptons: 

2 ^ — ^real + ^fake- 


The estimation of Vfake from MC simulation is difficult due to large uncertainties in the modeling of fake 
lepton production, and impractical due to a tiny fake rate requiring huge simulated samples. Therefore, 
the contribution of fakes was evaluated from data using the matrix method [48]. In addition to events 
from the signal data sample (labeled as “tight” events), a second (“loose”) set enriched with fake leptons 
was defined by removing the lepton isolation requirement. The number of events in each sample can be 
written as 

real take ’ 

^loose _ ^loose _|_ ^loose 
’^real ’^fake 


If the fractions of loose events that are also tight events are known, the number of events with fake leptons 
in the signal sample can be determined from a linear system of two equations. Given the probabilities 
for real and fake leptons that already passed the loose selection to also pass the tight selection, gi-eai = 
and gfake - ^feke'/^fakr’ number of loosc and tight events, the number of tight 

events with a fake lepton is 


tight _ 


■'frake 


Geal £fake 


£fake- 


The probability ffreai was determined with a tag-and-probe method [40] based on the identification of a 
tight lepton and a loose lepton in events originating from Z —> eejpp decays. The fake lepton probability 

^ When the selection is applied to electrons as well as muons, the thresholds refer to pj for muons and for electrons. 
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Cfake in the //+jets channel was derived for muons with large track do significance.'^ In the e+jets channel, 
Cfake was determined in the < 20 GeV, mxCVk) < 20 GeV region, where the number of events 

was corrected for the presence of real leptons. In both channels, the real and fake probabilities were 
parameterized in terms of lepton pj and 77 , the leading jet px. the distance between the lepton and the 
closest jet {AR(j), the ratio of the pj of the closest jet to the total number of jets, and the number 
of f 7 -tagged jets. The product of one-dimensional parameterizations was used, since there are too few 
events for a multidimensional parameterization. Based on this procedure, each event in the loose set was 
assigned a weight 

[(Geai - l)efake/(£reai “ £fake) for cvcnts with tight Icptons and ereaiSfake/CGeai “ ^fake) for the rest], and 
the weighted sample of data events selected using the loose object criteria was taken as the multijet 
contribution to the LHD background distribution. 

The numbers of observed and expected events are shown in Table 1. Uncertainties on the numbers of 
expected ti, single top, and diboson events are derived from theoretical uncertainties on the respective 
production cross-sections. The uncertainty on U-rjets events is evaluated based on a +4% uncertainty on 
inclusive U-i-jets production and a ±24% uncertainty on each additional jet [49], giving a total uncertainty 
of ±45% for U±jets background. This large uncertainty does not have an impact on the final result, 
since the normalization of the U±jets background is obtained from the data as described in Sec. 5. The 
uncertainty on the multijet background is evaluated by varying the parameters of the matrix method as 
described in Sec. 7. 


Table 1: Observed numbers of events in the e±jets and yu±jets channel together with estimated contributions from 
the tt signal and various background sources and associated uncertainties as described in the text. 


Event counts 

e±jets 

ju±jets 

tt 

1230001^^"^ 


U±jets 

38000±17000 

47000±21000 

Single top 

12100±600 

14900±800 

Dibosons 

710±40 

880±50 

Multijets 

2800±2400 

1500"^^°° 

Total 

^^PTTQqq+TOTD^ 

216000±23000 

Data 

176286 

220369 


5 Determination of the tt cross-section 

The number of tt events in each channel and their combination was obtained from the data using a template 
fit to the LHD distribution. The LHD function was constructed using the projective likelihood method 
defined in Ref. [50]. The LHD for event i is defined as the ratio of the signal L® to the sum of signal and 
background likelihood functions L® ± L|’ : 


Di = 


I 

i i 


do is the distance of closest approach in the transverse plane of a track to the primary vertex of the event; the do significance 
is the do value divided by its fitted uncertainty. 
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where the likelihood functions = n / = 11 / the products of the probability density func- 

tions of kinematic variables vj for signal and background, respectively. The variables used in the LHD 
were lepton pseudorapidity r]{ and transformed aplanarity Jl' = exp(- 8 j?l). The aplanarity ^ is defined as 
3/2 times the smallest eigenvalue of the momentum tensor pu^pj^j ^2^ where pik is 

the /-th momentum component and is the modulus of the momentum of object k. The objects included 
in the sum are leptons, jets, and The choice of variables followed the analysis in Ref. [51] and was 
motivated by optimization in terms of best separation between the signal and the background, the quality 
of the MC description of the variables in the data, and reduced sensitivity to the jet-energy scale. As the 
fraction of signal events in the selected data sample is large, two discriminating variables were found to 
be sufficient to provide an adequate signal-to-background separation. 

The and J{' distributions for data and MC simulation are shown in Fig. 1 . The correlation between these 
variables is found to be negligible in both signal and background. The LHD distributions for signal and 
background events are compared with the data in Fig. 2. The templates of LHD in the e-i-jets and p-tjets 
channels are different due to the fact that electrons in the barrel-endcap transition region are excluded and 
have lower efficiency in the forward region, whereas the efficiency as a function of 77 ^ is quite smooth. 

A weighted sum of the templates for the tt signal and all the backgrounds was fitted to a binned LHD 
distribution in the data. The contributions from single top and diboson events were fixed according to 
their theoretical cross-sections. The contribution from multijets events was obtained using the matrix 
method. 

The normalizations of the signal and F-i-jets (the dominant background) were left as free parameters in 
the fit. In each channel j, the fit was performed by minimizing the negative log-likelihood function for a 
Poisson model, 

L; ^ -In Lj = 2 J] (vij - riij + nij In ntj - nij In y,y), 
i 

summed over all bins i, where is the observed number of events and is the expected number of 
events in each bin or channel. The latter is defined as 

Vij = p“sij + p]bij + qtj, 

where sij, btj, and qtj are the predicted numbers of events for the tt signal, F-i-jets, and small backgrounds 
(single top, dibosons, and multijets), respectively, and p‘^ and pj are the parameters of the fit. The 
resulting value of the fit parameter p‘‘ was used to extract the number of tt events in the data. Figure 3 
shows the distribution of the LHD for data and MC simulation. In both Figs. 1 and 3, the contributions 
from tt and F-i-jets production are normalized according to the results of the LHD fit. 

The channels were combined by minimizing L = Yjj Lj where the sum is taken over the e-i-jets and p-i-jets 
channels. 

After the number of tt events in data Nti was obtained from the fit, the tt production cross-section in each 
channel was determined as 


o-tt 


Nti 


EtiXSx £’ 


( 1 ) 


where £tt is the signal-reconstruction efficiency (4.4% in the e-i-jets channel and 5.5% in the p-i-jets chan¬ 
nel), S is the tt —>> 1 lepton-i-jets branching ratio, and £ is the total integrated luminosity. 
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Figure 1; Distributions of the lepton pseudorapidity (top) and the modified aplanarity (bottom) in the e+jets (left) 
and /^+jets (right) channel. The data (dots) are compared to the SM expectation broken down into contributions from 
ff, single top, V+jets, diboson, and multijet production. The hatched areas correspond to the combined statistical 
and systematic uncertainties excluding the uncertainty on the ff signal modeling with different MC generators (see 
Sec. 7 for details). The bottom part of each plot shows the ratio of the data to the predicted value together with 
combined statistical and systematic uncertainties. 
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Figure 2: The templates of the likelihood discriminant LHD in the e+jets (left) and yU+jets (right) channels for 
signal, y+jets background, and fixed backgrounds, which include single top, diboson and multijets background. 
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Figure 3: Distributions of the likelihood discriminant LHD in the e+jets (left) and ju+jets (right) data and a weighted 
sum of templates from the ti signal and various backgrounds. The contributions from ti and y+jets production are 
normalized according to the results of the likelihood discriminant fit, the single top and dibosons contributions are 
normalized according to their theoretical cross-sections, and the multijet background is normalized using the data- 
driven method explained in the text. The bottom part of each plot shows the ratio of the data to the predicted value 
together with combined statistical and systematic uncertainties. 
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6 Fiducial cross-section measurement 


The tt cross-section in the fiducial volume was measured to allow for a more robust comparison to the 
theoretical prediction without extrapolating to regions outside of the detector acceptance. The definition 
of the fiducial volume is based on MC simulation and uses particle-level objects constructed using stable 
particles with a mean lifetime Tpartide > 0.3 x 10~^° s. Electrons and muons were required to originate 
from t Wb £vb decays, either directly or via a leptonically decaying r, and to have px > 40 GeV 
and \t]\ <2.5. The lepton momenta were corrected by adding the energy and momentum of photons inside 
a cone of radius AR = 0.1 around the lepton direction. Jets were created by clustering particles using the 
anti-A:f algorithm with radius parameter R = 0.4 (neutrinos, electrons and muons from IT-boson decays 
were excluded from particle-jet formation). Particles from the underlying event were included in this 
definition, but particles resulting from pileup were not. The particle-jets were required to have \t]\ < 2.5 
and pj > 25 GeV. The particle-jet f7-identification was based on the presence of nearby B-hadrons with 
a pt of at least 5 GeV using the ghost tagging method [52]. The overlap removal procedure included the 
removal of particle-jets within AR = 0.2 of the nearest electron with Ej > 25 GeV followed by removal 
of electrons and muons within AR = 0.4 of the nearest remaining particle-jet with px > 25 GeV. The 
events were required to have exactly one lepton with px > 40 GeV, no second lepton with px > 25 GeV, 
and at least three particle-jets, of which at least one jet was identified as a b-]&t. Missing transverse 
momentum was calculated using neutrinos from W leptonic decays. Similar to the reconstructed event 
selection, > 30 GeV and nix(lT) > 30 GeV requirements were applied. 

The simulated tf events were split into two categories: in-fiducial (satisfying the fiducial selection criteria) 
and out-of-fiducial (the rest). 

The number of tt events Nti obtained with the LHD fit was converted to the number of tt events in the 
fiducial volume using the fraction of reconstructed MC events passing the fiducial selection. The 
number of tt events in the fiducial volume was further used to determine the fiducial cross-section 
as 


fid 


CTp, 




( 2 ) 


where is the fiducial signal-reconstruction efficiency. The latter is defined as the ratio of the number 
of events passing both the reconstruction requirements and the particle-level selection, ^ the total 
number of events passing the particle level selection, N^'^\ 




Nl 


fid 




It was found to be 36% in the e-i-jets channel and 44% in the p-i-jets channel. 


7 Systematic uncertainties 

The systematic uncertainties in the tt cross-section measurement can be split into those related to the ob¬ 
ject reconstruction and momentum and energy measurements, those related to background evaluation, and 
those related to the signal modeling. The most important contribution to the total systematic uncertainty 
arises from effects related to the tt modeling, as discussed below. 
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To avoid statistical fluctuations, all systematic uncertainties were evaluated using ensemble tests. Simu¬ 
lated data sets (“ensembles”) were generated by picking random combinations of events from a pool of 
MC original (nominal) events according to the expected number of events due to signal and each back¬ 
ground type. For each source of systematic uncertainty, the modified events that passed the regular event 
selection were used to construct the LHD templates used in the ensemble tests. A weighted sum of mod¬ 
ified signal and background templates was fit to the LHD distributions for each ensemble. The results of 
the fits were averaged over ensembles using a Gaussian fit of the difference between the average number 
of tt events for the nominal and the modified template fits, and the fitted difference was propagated to the 
cross-section uncertainty. The MC statistical uncertainty associated to the limited signal and background 
sample size was found to be below ±0.1%. 

Lepton reconstruction: the uncertainties due to lepton trigger, identification, energy or momentum res¬ 
olution and reconstruction efficiencies were estimated from Z —> cc//iyu, Jjip eel jin, and W ^ ev 
processes using techniques discussed in Refs. [40,53,54]. These uncertainties are relatively small, dom¬ 
inated by the lepton identification in the c±jets channel (±2%) and the muon triggering efficiency in the 
/r±jets channel (±3%). They are specific to each lepton flavor and therefore uncorrelated between the 
channels. 

Jet reconstruction: the uncertainty on the tt cross-section due to the uncertainty on the jet-energy scale 
was estimated by varying the jet energies according to the uncertainties derived from simulation and in- 
situ calibration measurements using a model with 22 orthogonal components [44,55]. The variations in 
jet energies were also propagated to the value. The uncertainty due to the difference in jet-energy 
resolution between the data and MC events was evaluated by smearing the MC jet transverse momentum 
according to the jet resolution as a function of the jet pT and p [56]. The uncertainty due to the jet 
reconstruction efficiency was estimated by randomly discarding jets according to the difference in jet 
reconstruction efficiency between the data and MC. The lower value of jet vertex fraction was varied 
between 0.4 and 0.6 as motivated by the Z —> cc/pp±jets studies [57]. The dominant jet reconstruction 
systematic uncertainty on the tt cross-section is due to the uncertainty on the jet-energy scale (±3%). 

^-tagging: the uncertainties on the (7-tagging scale factors were obtained separately for (7-jets, c-jets and 
light flavor jets, and were treated as uncorrelated for each type of jet. The uncertainties on the inefficiency 
scale factors, applied when a jet was not tagged, were treated as fully anticorrelated to those from the 
efficiency scale factors. Since the efficiency to tag a (7-jet is 70%, and only one of the two ( 7 -jets in the 
event is required to be tagged, the tt event tagging efficiency is large (about 90%) and its uncertainty 
is dominated by ( 7 -jet tagging uncertainties, which are small compared to non-(7-jet ones. The resulting 
( 7 -tagging contribution to the overall systematic uncertainty is about ±2%. 

Missing transverse momentum: the systematic uncertainties associated with the momenta and energies 
of reconstructed objects (leptons and jets) were also propagated to the calculation. The recon¬ 
struction also receives contributions from the presence of low-px jets and calorimeter cells not included 
in reconstructed objects (“soft terms”). The systematic uncertainty on soft terms was evaluated using 
Z ^ pp events from the data/MC ratio in events without jets and from the balance between soft 
terms and hard objects using methods similar to those used in Ref. [58]. The measurement is not a 
significant source of systematic uncertainty (below ±1%). 

W+jets model: this source of systematic uncertainty does not affect the number of tt events expected after 
applying the selection requirements but it does change the shape of the background template. This uncer¬ 
tainty was estimated using templates with lT±jets events generated by Sherpa instead of Alpgen±Pythia 


12 


and was found to be small (below 1%) and to partially cancel out when combining the e+jets and /i+jets 
channels. 

W + heavy-flavor composition: the heavy-flavor composition of W-i-jets events was studied using W 
events with two additional jets split into and samples.^ Uncertainties were derived from contribu¬ 
tions from W+c,W +CC, and W+bb based on the statistical uncertainty, uncertainties on the extrapolation 
from 2-jets to >3-jets, and MC modeling, following the procedure described in Ref. [60]. The uncertain¬ 
ties themselves are very large (about ±50%), but the LHD template shapes for IT-bosons produced in 
association with different jet flavors are similar and since the normalization of the lT±jets background 
is derived from a fit to the data, the resulting systematic uncertainty on the tt production cross-section is 
small (±1%). 

Small backgrounds: uncertainties on single top-quark and diboson production were evaluated by varying 
the relevant theoretical cross-sections within their uncertainties ±lo-theor. which resulted in a tt produc¬ 
tion cross-section uncertainty of the order of ±1%. The uncertainty on the multijet contribution was 
derived from a comparison with results with fireai obtained using data control samples with large mj{W) 
(ju±jets) or large (e±jets) depleted of fake leptons; the variation of gfake due to the change of the low 
region; and the choice of MC samples used to subtract the real leptons. As seen in Table 1, 
the uncertainty on the number of multijet events is very large (about ±100%), mainly due to uncertainties 
in efake^ but their fraction in the data sample is small due to the high lepton-pj requirement, so the overall 
effect is about ±1%. 

Pileup: the MC samples used in the analysis were reweighted to reproduce correctly the distribution of 
the number of pp collisions per bunch crossing in the data. The accuracy of the average correction factor 
due to pileup reweighting was found to be ±4%. The pileup reweighting uncertainties originated from 
the modeling of pileup events, including uncertainties in the pp inelastic cross-section. The resulting 
effect on the measured tt cross-section is small (about ±0.1%) but is still included in the total systematic 
uncertainty. 

Initial- and flnal-state radiation (ISR and FSR) modeling: ISR or FSR changes the number of jets 
in the event. To evaluate the uncertainty linked to the modeling of the ISR or FSR, tt MC samples with 
modified ISR and FSR modeling were used. The MC samples used for the evaluation of this uncertainty 
were generated using the AcerMC generator interfaced to Pythia, where the parameters of the generation 
(Aqcd> Gmax scale, transverse momentum scale for space-like parton-shower evolution [20]) were varied 
to span the ranges compatible with the results of a measurement of tt production with a veto on additional 
central jet activity [61]. This uncertainty is large for the total inclusive tt production cross-section (±3%) 
but significantly reduced (below ±1%) for the fiducial cross-section measurement. 

MC generator: the choice of MC generator used in the signal modeling affects the kinematic properties 
of simulated tt events and reconstruction efficiencies. For the purpose of addressing this effect, tt events 
simulated with MC@NLO interfaced to Herwig±Jimmy were used and compared to results obtained with 
PowHEG interfaced to Herwig±Jimmy. The resulting systematic uncertainty was found to be ±3%. 

MC parton shower and fragmentation model: the effect of the parton shower and fragmentation mod¬ 
eling was estimated by comparing the results obtained from the default MC sample simulated by Powheg 
interfaced to Pythia with a modified model simulated by Powheg interfaced to Herwig±Jimmy. The 


® The ratio of W* to W production cross-sections at the LHC is predicted much more precisely than the cross-sections them¬ 
selves [59]. 
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change in the number of expected tt events due to the modified selection efficiency and the LHD dis¬ 
tribution was propagated to the systematic uncertainty on the tt production cross-section. Using this 
approach, the systematic uncertainty due to the MC parton shower and fragmentation model was found 
to be ±2%. 

MC parton distribution functions (PDFs): this uncertainty is one of the most significant ones in this 
analysis (±6%). An event-by-event reweighting was applied to the MC@NLO tt sample generated us¬ 
ing the central value of the CTIO PDF. The CTIO variations, as well as both the central values and the 
variations of MSTW2008 NLO at the 68% confidence level and NNPDF23 were used fo evaluafe fhe 
sysfemafic uncerfainfy, following fhe PDF4LHC recommendations [9]. For each PDF and ifs variafions 
fhe change in fhe tt cross-secfion was calculafed using ensemble fesfs. The final envelope of all variafions 
was symmefrized and half of ifs size was quoted as fhe sysfemafic uncerfainfy due fo fhe PDF uncer¬ 
fainfy. The effecf arises primarily from fhe gluon PDF componenf, which is large for signal and small for 
background. Therefore fhe PDF uncertainly is dominated by fhe choice of PDF in fhe tt simulation. The 
PDF uncerfainfy affecls bolh fhe LHD shape and fhe selection efficiency, and is significanl for bofh fhe 
inclusive and fhe fiducial cross-secfion measuremenl. 

Luminosity: this uncertainty was evaluated separately and not combined with the rest of the systematic 
uncertainties. Following the method described in Ref. [62], the uncertainty on the total integrated lumi¬ 
nosity was determined to be ±2.8% from a preliminary calibration of the luminosity scale derived from 
beam-separation scans performed in November 2012. By combining the change in X in Eqs. (1) and (2) 
and the change in the contribution of single top and diboson events, the overall uncertainty due to the 
integrated luminosity was found to be ±2.9%. 

LHC beam energy: calibrations determined the beam energy to be (0.30 ± 0.66)% smaller than the 
nominal value of 4 TeV per beam [63]. The measured tt cross-section value is not corrected for this shift. 
However, an uncertainty of ±1.7%, corresponding to the expected change in cTtf for a ±0.66% change 
in ^/s is quoted separately in the final result. The effect of changing the LHC beam energy on Stf is 
negligible. 

The evaluation of the systematic uncertainties other than the rt MC dependence for the measurement in the 
fiducial volume is the same as for the inclusive measurement. For the systematic uncertainties obtained 
using different tt MC samples, the different signal-reconstruction efficiency in the fiducial volume was 
taken into account. 

The systematic uncertainties due to the jet reconstruction and measurement are grouped together, 
as well as the systematic uncertainties due to the background evaluation. A summary of the systematic 
uncertainties is given in Table 2. 

8 Results and conclusions 

Inclusive and fiducial ti production cross-sections have been measured at the LHC in pp collisions at 
- 8 TeV using data collected with the ATLAS detector in 2012, corresponding to 20.3 fb“^ of inte¬ 
grated luminosity. The measurements were in the lepton±jets final state using a likelihood discriminant 
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Table 2; Summary of the systematic uncertainties in the measurements of the ff production cross-section (%). The 
bottom part of the table lists systematic uncertainties for the fiducial ff production cross-section for the cases where 
they are different from those for the inclusive one. 


Uncertainty on inclusive <Ttt 

e-i-jets 

/r+jets 

t?+jets 

Lepton reconstruction 

+2.7 

-2.6 

+2.1 

-1.9 

+ 1.7 

-1.6 

Jet reconstruction and 

+3.3 

-3.9 

+2.6 

-3.2 

+2.8 

-3.4 

fi-tagging 

+2.1 

-1.9 

+2.2 

-1.9 

+2.1 

-1.9 

Backgrounds 

+2.8 

-3.0 

+ 1.8 

-2.1 

+ 1.7 

-2.1 

Monte Carlo generator 

-2.2 

+2.2 

-3.3 

+3.3 

-2.7 

+2.7 

Parton shower and fragmentation 

+2.0 

-2.0 

+2.6 

-2.6 

+2.3 

-2.3 

Initial- and final-state radiation 

-4.1 

+4.1 

-1.8 

+1.8 

-3.0 

+3.0 

Parton distribution functions 

+6.2 

-6.0 

+5.6 

-5.9 

+5.9 

-5.9 

Total 

+9.7 

-9.8 

+8.4 

-8.7 

+8.6 

-8.9 

Uncertainty on fiducial cr„- 

e+jets 

ju+jets 

T’+jets 

Monte Carlo generator 

-2.1 

+2.1 

-3.5 

+3.5 

-2.8 

-2.8 

Parton shower and fragmentation 

-2.6 

+2.6 

-3.1 

+3.1 

-2.9 

+2.9 

Initial- and final-state radiation 

+0.4 

-0.4 

+0.2 

-0.2 

+0.3 

-0.3 

Total 

+8.9 

-9.0 

+8.5 

-8.8 

+8.3 

-8.6 


fit. Assuming a top-quark mass of 172.5 GeV, the inclusive tt production cross-section is found to be 

e-i-jets : (Tti = 256 ± 2(stat.) ± 25(syst.) ± 7(lumi.) ± 4(beam) pb, 

//-i-jets : (Ttj = 260 ± l(stat.) - 8(lumi.) ± 4(beam) pb, 

T’-f-jets : (Ttt = 258 ± l(stat.) il^^sCsyst.) ± 8(lumi.) ± 4(beam) pb, 

where the four quoted uncertainties are due to the number of events in data (statistical uncertainty), 
systematic effects (as described in Sec. 7), the imprecisely known integrated luminosity, and the LHC 
beam energy. Since the results were obtained using a fixed value of mtop, a set of ti samples with different 
top-quark masses was generated in order to study the mtop dependence of the measured tt production 
cross-section, which is found to be (Acr,f/crff)/Amtop = -1.1 %/GeV in all channels. 

The fiducial tt cross-section is found to be 

e-i-jets : = 11.3 ± 0.1(stat.) ± 1.0(syst.) ± 0.3(lumi.) ± 0.2(beam) pb, 

/r-i-jets : = 11.5 ± 0.1(stat.) ± 1.0(syst.) ± 0.3(lumi.) ± 0.2(beam) pb, 

T’-i-jets : = 22.8 ± 0.1(stat.) - 0.7(lumi.) ± 0.4(beam) pb, 

where the four quoted uncertainties represent the same sources as indicated before. 

The largest systematic uncertainties are due to ff MC modeling, including the PDF uncertainty, the choice 
of MC generator, the parton shower and fragmentation model and initial- and final-state radiation. The 
systematic uncertainty from ISR or FSR is found to be reduced for the fiducial cross-section measurement 
compared to the inclusive case, while the uncertainties due to the choice of MC generator and the parton 
shower and fragmentation model are found to increase slightly in the fiducial cross-section measurement 
compared to the inclusive case. 
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The measured inclusive ti production cross-section is in good agreement with the NNLO-i-NNLL theory 
prediction, cr,,- = 253^}^ pb [2]. It is also consistent with the ATLAS and CMS measurements of the 
tt production cross-section in the dilepton channel [15,16]. The ATLAS measurement had the result of 
cr,f = 242.4 ± L7(stat.) ±5.5(syst.) ±7.5(lumi.) ± 4.(beam) pb [15], where the four quoted uncertainties 
were due to the same sources as above. 
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